Fast Converging Path Integrals for Time-Dependent Potentials I: Recursive Calculation of Short-Time Expansion of the Propagator

In this and subsequent paper arXiv:1011.5185 we develop a recursive approach for calculating the short-time expansion of the propagator for a general quantum system in a time-dependent potential to orders that have not yet been accessible before. To this end the propagator is expressed in terms of a discretized effective potential, for which we derive and analytically solve a set of efficient recursion relations. Such a discretized effective potential can be used to substantially speed up numerical Monte Carlo simulations for path integrals, or to set up various analytic approximation techniques to study properties of quantum systems in time-dependent potentials. The analytically derived results are numerically verified by treating several simple models.Comment: 29 pages, 5 figure

Faraday waves in binary non-miscible Bose-Einstein condensates

We show by extensive numerical simulations and analytical variational calculations that elongated binary non-miscible Bose-Einstein condensates subject to periodic modulations of the radial confinement exhibit a Faraday instability similar to that seen in one-component condensates. Considering the hyperfine states of $^{87}$Rb condensates, we show that there are two experimentally relevant stationary state configurations: the one in which the components form a dark-bright symbiotic pair (the ground state of the system), and the one in which the components are segregated (first excited state). For each of these two configurations, we show numerically that far from resonances the Faraday waves excited in the two components are of similar periods, emerge simultaneously, and do not impact the dynamics of the bulk of the condensate. We derive analytically the period of the Faraday waves using a variational treatment of the coupled Gross-Pitaevskii equations combined with a Mathieu-type analysis for the selection mechanism of the excited waves. Finally, we show that for a modulation frequency close to twice that of the radial trapping, the emergent surface waves fade out in favor of a forceful collective mode that turns the two condensate components miscible.Comment: 13 pages, 10 figure

Development of Grid e-Infrastructure in South-Eastern Europe

Over the period of 6 years and three phases, the SEE-GRID programme has established a strong regional human network in the area of distributed scientific computing and has set up a powerful regional Grid infrastructure. It attracted a number of user communities and applications from diverse fields from countries throughout the South-Eastern Europe. From the infrastructure point view, the first project phase has established a pilot Grid infrastructure with more than 20 resource centers in 11 countries. During the subsequent two phases of the project, the infrastructure has grown to currently 55 resource centers with more than 6600 CPUs and 750 TBs of disk storage, distributed in 16 participating countries. Inclusion of new resource centers to the existing infrastructure, as well as a support to new user communities, has demanded setup of regionally distributed core services, development of new monitoring and operational tools, and close collaboration of all partner institution in managing such a complex infrastructure. In this paper we give an overview of the development and current status of SEE-GRID regional infrastructure and describe its transition to the NGI-based Grid model in EGI, with the strong SEE regional collaboration.Comment: 22 pages, 12 figures, 4 table

Stability of Intercelular Exchange of Biochemical Substances Affected by Variability of Environmental Parameters

Communication between cells is realized by exchange of biochemical substances. Due to internal organization of living systems and variability of external parameters, the exchange is heavily influenced by perturbations of various parameters at almost all stages of the process. Since communication is one of essential processes for functioning of living systems it is of interest to investigate conditions for its stability. Using previously developed simplified model of bacterial communication in a form of coupled difference logistic equations we investigate stability of exchange of signaling molecules under variability of internal and external parameters.Comment: 11 pages, 3 figure

Geometric Resonances in Bose-Einstein Condensates with Two- and Three-Body Interactions

We investigate geometric resonances in Bose-Einstein condensates by solving the underlying time-dependent Gross-Pitaevskii equation for systems with two- and three-body interactions in an axially-symmetric harmonic trap. To this end, we use a recently developed analytical method [Phys. Rev. A 84, 013618 (2011)], based on both a perturbative expansion and a Poincar\'e-Lindstedt analysis of a Gaussian variational approach, as well as a detailed numerical study of a set of ordinary differential equations for variational parameters. By changing the anisotropy of the confining potential, we numerically observe and analytically describe strong nonlinear effects: shifts in the frequencies and mode coupling of collective modes, as well as resonances. Furthermore, we discuss in detail the stability of a Bose-Einstein condensate in the presence of an attractive two-body interaction and a repulsive three-body interaction. In particular, we show that a small repulsive three-body interaction is able to significantly extend the stability region of the condensate.Comment: 27 pages, 13 figure

Computation of Asteroid Proper Elements on the Grid

A procedure of gridification of the computation of asteroid proper orbital elements is described. The need to speed up the time consuming computations and make them more efficient is justified by the large increase of observational data expected from the next generation all sky surveys. We give the basic notion of proper elements and of the contemporary theories and methods used to compute them for different populations of objects. Proper elements for nearly 70,000 asteroids are derived since the beginning of use of the Grid infrastructure for the purpose. The average time for the catalogs update is significantly shortened with respect to the time needed with stand-alone workstations. We also present basics of the Grid computing, the concepts of Grid middleware and its Workload management system. The practical steps we undertook to efficiently gridify our application are described in full detail. We present the results of a comprehensive testing of the performance of different Grid sites, and offer some practical conclusions based on the benchmark results and on our experience. Finally, we propose some possibilities for the future work

Dynamics of weakly interacting bosons in optical lattices with flux

Realization of strong synthetic magnetic fields in driven optical lattices has enabled implementation of topological bands in cold-atom setups. A milestone has been reached by a recent measurement of a finite Chern number based on the dynamics of incoherent bosonic atoms. The measurements of the quantum Hall effect in semiconductors are related to the Chern-number measurement in a cold-atom setup; however, the design and complexity of the two types of measurements are quite different. Motivated by these recent developments, we investigate the dynamics of weakly interacting incoherent bosons in a two-dimensional driven optical lattice exposed to an external force, which provides a direct probe of the Chern number. We consider a realistic driving protocol in the regime of high driving frequency and focus on the role of weak repulsive interactions. We find that interactions lead to the redistribution of atoms over topological bands both through the conversion of interaction energy into kinetic energy during the expansion of the atomic cloud and due to an additional heating. Remarkably, we observe that the moderate atomic repulsion facilitates the measurement by flattening the distribution of atoms in the quasimomentum space. Our results also show that weak interactions can suppress the contribution of some higher-order nontopological terms in favor of the topological part of the effective model

Consequences of increased longevity for wealth, fertility, and population growth

We present, solve and numerically simulate a simple model that describes the consequences of increased longevity on fertility rates, population growth and the distribution of wealth in developed societies. We look at the consequences of the repeated use of life extension techniques and show that they represent a novel commodity whose introduction will profoundly influence key aspects of economy and society in general. In particular, we uncover two phases within our simplified model, labeled as 'mortal' and 'immortal'. Within the life extension scenario it is possible to have sustainable economic growth in a population of stable size, as a result of dynamical equilibrium between the two phases.Comment: 13 pages, 5 figures, uses elsart.cl

Dyons in Nonabelian Born-Infeld Theory

We analyze a nonabelian extension of Born--Infeld action for the SU(2) group. In the class of spherically symmetric solutions we find that, besides the Gal'tsov--Kerner glueballs, only the analytic dyons have finite energy. The presented analytic and numerical investigation excludes the existence of pure magnetic monopoles of 't Hooft--Polyakov type.Comment: 12 pages, 3 figure

Fast Convergence of Path Integrals for Many-body Systems

We generalize a recently developed method for accelerated Monte Carlo calculation of path integrals to the physically relevant case of generic many-body systems. This is done by developing an analytic procedure for constructing a hierarchy of effective actions leading to improvements in convergence of $N$-fold discretized many-body path integral expressions from 1/N to $1/N^p$ for generic $p$. In this paper we present explicit solutions within this hierarchy up to level $p=5$. Using this we calculate the low lying energy levels of a two particle model with quartic interactions for several values of coupling and demonstrate agreement with analytical results governing the increase in efficiency of the new method. The applicability of the developed scheme is further extended to the calculation of energy expectation values through the construction of associated energy estimators exhibiting the same speedup in convergence.Comment: 20 pages, 6 figures, 1 tabl